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ABSTRACT 

We present low- and high-resolution mid-infrared spectra and photometry for eight compact sym- 
metric objects (CSOs) taken with the Infrared Spectrograph on the Spitzer Space Telescope. The 
hosts of these young, powerful radio galaxies show significant diversity in their mid-IR spectra. This 
includes multiple atomic fine-structure lines, H2 gas, PAH emission, warm dust from T = 50 to 150 K, 
and silicate features in both emission and absorption. There is no evidence in the mid-IR of a single 
template for CSO hosts, but 5/8 galaxies show similar moderate levels of star formation (< IOMq yr _1 
from PAH emission) and silicate dust in a clumpy torus. The total amount of extinction ranges from 
Ay ~ 10 to 30, and the high-ionization [Ne V] 14.3 and 24.3 /zm transitions are not detected for 
any galaxy in the sample. Almost all CSOs show contributions both from star formation and active 
galactic nuclei (AGNs), suggesting that they occupy a continuum between pure starbursts and AGNs. 
This is consistent with the hypothesis that radio galaxies are created following a galactic merger; 
the timing of the radio activity onset means that contributions to the infrared luminosity from both 
merger-induced star formation and the central AGN are likely. Bondi accretion is capable of powering 
the radio jets for almost all CSOs in the sample; the lack of [Ne V] emission suggests an advection- 
dominated accretion flow (ADAF) mode as a possible candidate. Merging black holes (BHs) with 
Mbh > 10 8 Mq likely exist in all of the CSOs in the sample; however, there is no direct evidence from 
this data that BH spin energy is being tapped as an alternative mode for powering the radio jets. 
Subject headings: accretion - accretion disks - galaxies: evolution - galaxies: jets - infrared: galaxies 
- radio continuum:galaxies 



1. INTRODUCTION 

An enduring mystery in the study of radio galax- 
ies is the nature of their power source. Specifi- 
cally, most galaxies which possess a radio-luminous 
active galactic nucleus (AGN) and rapidly advancing 
jets show little evidence for a central accreting disk. 
For example, in Centaurus A and Virgo A (M87) 
there is no blue continuum or broad emission lines, 
but also no IR-emitting screen large enough to ob- 
scure a direct view of a typical broad emission line re- 
gion dPerlman et aL 2001a; Whvsong fc Antonuccill2004l ; 
iRadomski et al.ll2008l ). While this modern paradigm uni- 
fies all AGN as supermassive black holes (BH) fed by 
an accretion disk orbiting the BH, many AGN show lit- 
tle evidence for the large amounts of circumnuclear gas 
needed to power them. These and other considerations 
led to the proposal that some AGN could be dominated 
by "advection-domi nated accretion Sows" (ADAFs; e.g. 
INaravan fc Yll[l99l . 

An alternative is that the active phase may be trig- 
gered by a major merger of disk galaxies and, con- 
currently, a merger of their supermassive BHs. It is 
now understood on the basis of detailed numerical sim- 
ulations that the merger of disk g alaxies can create a 
relati vely normal elliptical galaxy ()Barnes fc Hernquistl 
1992). Since most (if not all) bright disk galaxies contain 
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supermassive BHs, Wil son fc Colbert! (|1995l) proposed a 
scenario for the creation of a radio-loud AGN in which 
the merger of two disk galaxies and their BHs can spin up 
the resulting BH, thus powering the AGN and its jets. 
The detection of X-ray point sources in radio-loud el- 
lipticals provi des observational evidence for the mergin g 
BH scenario (|Komossa et al.ll2003l : iBianchi et al.| [2008). 
BH spin accounts naturally for the diminished output of 
successive radio outbursts in the brightest cluster galaxy 
(BCG) because the spin energy is diminished by each 
outburst and the accretion rate is too low for replenish- 
ment. At high redshift a merger of disk galaxies (a so- 
called "wet merger") and their individual supermassive 
BHs creates a rapidly-spinning supermassive BH and a 
luminous quasar. As the cluster or group of galaxies in 
which the BCG is imbedded develops a hot, dense ICM, a 
large fraction of cluster galaxies are stripped of their cold 
and warm gas. As a result, all such mergers are "dry" 
and provide no additional accretion power. The BH spin 
hypothesis thus predicts the correct sign of radio-loud 
AGN evolution, which is not explained by the simplest 
accretion scenarios. 

An ideal class of galaxies in which to test these compet- 
ing hypotheses are compact symmetric objects (CSOs). 
CSOs are radio-loud galaxies characterized by jet and/or 
hot spot activity on both sides of a central engine. They 
arc morphologically similar to classi cal double radio 
sources, but with size scales o f < 1 kpc (jPhillips fc Mute] 
119821 : iWilkinson et al.lll994f ). The radio emission typi- 
cally peaks near frequencies of a few GHz, which can be 
explained in mos t CSOs by syn c hrotro n self- absorption 
in the galaxies (|de Vries et al.l 12009 a. b). Kinematic 
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ages of the radio jets show they are extremely young , 
with ages of £ 40 00 yr (|Owsianik fc Conwavl H991 
IGughucci et al.ll2005ft . As a result, it has been suggested 
that CSOs represent the early st ages of an AGN evoly - 
ing into an FR II radio galaxy (Rcadhc ad et al.| [l996'), 
akin to a "mini-Cygnus A" (Bcgclnianl ll996D . Given the 
young ages of these galaxies, they offer a much better 
chance for detecting signs of a recent merger that may 
be powering the radio activity. 

HST imaging of t hree nearby CSO host galaxies 
(jPerlman et al.ll2001bft show them to be nearly normal el- 
lipticals; however, disturbed isophotes at the outer edges 
suggest that a past me rger has taken place in all objects 
(|Perlman et al.l l2~001bD . Simulations suggest that the 
radius of the photometric irregularities determines the 
timescale since t he merger due to virilization of the stel- 
lar d istribution (jMihos et al.l 119951: iMihos fe Hernquistl 
1996). This is evidence for a recent merger in all 
three CSOs which commenced ~ 100 million years ago, 
the time delay require d for supermassive BHs to merge 
(|Begelman et al.l fl984f) . This is also comparable to the 
timescale for driving gas to the center of the merger sug- 
gested by other simulations (iDi Matteo et alT 2005. 2008: 
iHopkins et alll2005t iSpringel et al.ll2005aU bl). While this 
timescale supports the Wilson-Colbert hypothesis, it of- 
fers no concrete proof that radio-loud AGN are powered 
by BH spin. The accretion model requires large amounts 
of cold or warm gas close to the CSO, however, which is 
not seen - therefore, the power source for the jets is still 
an open question. 

Since distinguishing between a model for radio-loud 
AGN fueled by BH spin or one fueled by accretion will 
significantly affect our understanding of jet production 
and the differences between radio-loud and radio-quiet 
AGN, CSOs provide a unique opportunity to investigate 
the nuclear conditions in radio galaxies recently triggered 
by some mechanism into activity. One key question is 
whether large amounts of nuclear obscuration unresolved 
in the HST optical images could hide the accreting gas. 
Mid-IR observations should detect this gas if it exists, 
which may provide strong supporting evidence for ac- 
cretion flows leading to AGN outbursts in CSOs. Ob- 
servations like the HST imaging and the mid-IR results 
presented in this paper are the first that address such a 
hypothesis. 

To this end we observed a sample of nearby CSOs with 
the Spitzer Space Telescope with the goal of distinguish- 
ing between these two hypotheses and characterizing the 
properties of CSO host galaxies, including gas and dust 
content and their star formation rates. If accretion is 
the power source for these galaxies, then mid-IR pho- 
tons should penetrate any extinction screen that might 
obscure the AGN at shorter wavelengths; accreting gas 
should also emit via high-ionization atomic transitions 
due to the proximity of the AGN. Absence of this would 
suggest a "naked" AGN with no luminous accreting gas, 
supporting a model in which radio-loud AGN are pow- 
ered by BH spin due to either rapidly spinning BHs or 
binary BHs created in a merger. 

In |j2j we describe the CSO sample and the observa- 
tions taken with Spitzer, we describe the data reduction 
process in SQ] presents results of the observations, in- 
cluding full spectra and data tables for all measured fea- 
tures. We derive physical properties for the CSO hosts 



in iJ5] and classify them in the mid-IR. In [JS] we discuss 
the evidence for the power source of CSOs and their role 
in the evolution of radio galaxies. An Appendix summa- 
rizes the jet properties and multi-wavelength observa- 
tions of the host galaxies from the literature. We assume 
the WMAP5 cosmology with H = 70.5 km s" 1 Mpc -1 , 
Q M = 0.27, and Q A = 0.73. 

2. CSO SAMPLE 

Our goal for the sample was to select all known 
CSOs with z < 0.1. This was motivated by the need 
for the galaxies to be IR-bright (detectable in reason- 
able integration times with the IRS) and having large 
amounts of multi-wavelength data to characterize the 
hosts. Radio surveys have detected CSOs out to red- 
shifts of z ~ 0.5, typically confirming their CSO sta- 
tus t hrough VLBI (jPeck fe Tavlorl 12000 : iTavlor fe Peck 
2003 ) . In addition to objects selected at z < 0.1 JO'Dea 



1998; Pih lstrom et~aH 12001 iPolatidis fe Colrwavl l2003) 



we added the unusually IR-bright CSO PKS 1413+135, 
located at a redshift of z = 0.2467. Recent IDs of CSOs 
at z < 0.1 mean that the sample is likely not complete 
(e.g- IAugusto et al.ll2006ft . 

The galaxies were observed in two separate programs 
during Spitzer cycles 3 and 5. The CSO 4C +12 50 was 
previously observed with the IRS (|Farrah et al.1 12007ft 
and we used data publicly available through the Spitzer 
archive for this object. The original sample also included 
IRS observations of the Seyfert 2 galaxy NGC 5793. Al- 
though the gala xy does possess parse c-scale symmetric 
radio emission (|Hagiwara et al.l 12000ft , the radio spec- 
trum is not gigahertz-peaked and it has no proper mo- 
tion measurements that constrain the age. Data for this 
galaxy is thus not included in this paper. 

Table [T] lists the basic properties of the nine CSOs in 
the final sample. The heliocentric redshift is from opti- 
cal measurements with the exception of PK S 1413+135, 
which uses the redshift from Hi absorption (Cari flTet al.l 
1992). Measurements of the Hi column depend on the 
assumed spin temperature T s , typically on the order of 
~ 100 K in these galaxies. Optical classifications of the 
hosts are: weak- line radio galaxy (WLRG), narrow-line 
radio galaxy (NLRG), low- ionization nuclear emission 
line region (LINER), BL Lac, or broad-line radio galaxy 
(BLRG). The radio lumi nosity is the sp ecific power at 
1.4 GHz from the NVSS ([Condon et al.l fl998). The in- 
frared luminosity is from IRAS fluxes and the prescrip- 
tion of lSanders fc Mirabell (fl996l ). The X -ray luminosity 
is measured in the 2-10 keV band (see Appendix for ref- 
erences). 

3. SPITZER DATA REDUCTION 

The IRS contains six different modules in both low- 
(LR) and high-resolution (HR) modes (jHouck et al.1 
l200l . The short-high (SH) and long-high (LH) mod- 
ules operate at a resolution of R ~ 600 with wavelength 
coverage from 9.9 — 37.2 /im; the short-low (SL1 & SL2) 
and long-low (LL1 & LL2) modules operate at resolutions 
R ~ 56— 127, with coverage spanning 5.2 — 38 /im. There 
also exist two "bonus segments" that cover the overlap- 
ping range between the 1 st and 2 nd orders in both SL and 
LL; we used these to verify the flux calibration between 
the first and second orders. 

In addition to the spectroscopic modules, the IRS also 
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TABLE 1 

Properties of CSOs observed with Spitzer IRS 



Object 


RA 


Dec 


z s 


D L 


Hi column 


HI 


Optical 


Optical 


P ' y 0,(0,0 


log L IR 


log Lx 




J2000.0 


J2000.0 




[Mpc] 


[cm" 2 ] 


ref. 


spec. 


ref. 


[W/Hz] 


[L©] 


[L©] 


4C +31.04 


01 19 35.0 


+32 10 50 


0.0602 


264 


1.08 x 10 21 


(1) 


WLRG 


(8) 


25.34 


10.60 




4C +37.11 


04 05 49.2 


+38 03 32 


0.055 


242 


1.8 x 10 20 


(2) 


NLRG 


(9) 


25.05 






1146+59 


11 48 50.3 


+59 24 56 


0.0108 


48.2 


1.82 x 10 21 


(1) 


LINER 


(10) 


23.10 


9.09 




1245+676 a 


12 47 33.3 


+67 23 16 


0.1073 


495 


6.73 x 10 20 


(3) 


WLRG 


(8) 


25.02 






4C +12.50 


13 47 33.3 


+12 17 24 


0.1217 


571 


T s 6.2 x 10 18 


(4) 


NLRG 


(11) 


26.30 


12.15 


43.3 


OQ 208 


14 07 00.4 


+28 27 15 


0.0766 


349 


1.83 x 10 20 


(5) 


BLRG 


(12) 


25.08 


11.33 


42.7 


PKS 1413+135 


14 15 58.8 


+13 20 24 


0.2467 


1244 


4.6 x 10 22 


(6) 


BL Lac 


(13) 


26.31 


11.97 


44.4 


PKS 1718-649 


17 23 41.0 


-65 00 37 


0.0142 


62.6 




LINER 


(14) 


24.25 






1946+70 


19 45 53.5 


+70 55 49 


0.1008 


460 


2.2 x 10 23 


(7) 




(15) 


25.39 







References. (1) - Ivan Gorkom et all 119891) : (21 - IManess et all 120041) : (31 - ISaikia et all 120071) : (4) - [Mirabell fT98Sl); (5) - 

Vcrmculcn ct al (2003): (6) - Pcrlman ct al. (2002); (7) - Peck ct al (1999): (8) - Marcha ct al. (1996): (9) - Rodriguez et al. (2006); (10) 
-IKiml U989D ; (11 1 - iBaan et all 1199811 ; (121 - IPetrosian et all (120071) : ( 131 - Istocke et all (119921) ; (1 41 - iFilippenkol 11985ll ; (15) - [Snellen et all 
(119991) 

a Undetected with the spectral modules on the IRS; sec E|3.1I 



TABLE 2 
IRS observation log 



Object 


Date 


SL1 


SL2 


LL1 


LL2 


SH 


LH 


Program 


4C +31.04 


2007 Sep 01 


60 x 3 


60 x 3 


120 x 5 


120 x 5 


120 x 6 


240 x 4 


30515 


4C +37.11 


2009 Apr 10 


60 x 3 


60 x 3 


120 x 5 


120 x 5 


120 x 9 


240 x 11 


50591 


1146+59 


2007 May 01 


60 x 2 


60 x 2 


30 x 3 


30 x 3 


30 x 6 


60 x 2 


30515 


1245+676 


2009 Apr 19 


60 x 4 


60 x 4 


120 x 6 


120 x 6 


120 x 5 


240 x 7 


50591 


4C +12.50 


2004 Jan 07 


14 x 3 


14 x 3 


30 x 2 


30 x 2 


30 x 6 


60 x 4 


105 


OQ 208 


2007 Jun 24 


60 x 2 


60 x 2 


30 x 3 


30 x 3 


30 x 6 


60 x 3 


30515 


PKS 1413+135 


2007 Jul 31 


60 x 2 


60 x 2 


30 x 4 


30 x 4 


120 x 5 


60 x 6 


30515 


PKS 1718-649 


2009 Apr 18 


14 x 3 


14 x 3 


30 x 5 


30 x 5 


30 x 5 


60 x 5 


50591 


1946+70 


2009 Apr 07 


60 x 3 


60 x 3 


120 x 5 


120 x 5 


120 x 9 


240 x 12 


50591 



Note. — Exposure times for all modules are given as seconds per cycle X number of cycles. 



contains two peakup imaging arrays used to position the 
target on the spectral slit. We took dedicated 16 and 
22 /zm sample-up-the-ramp (SUR) peakups for each ob- 
ject, giving two photometric data points used to calibrate 
the absolute flux scale of the spectra. The dedicated 
peakup observations were not coincident in time with the 
spectroscopy. Half of our sources are unresolved in the 
IRS slit and can be treated as point sources; for the re- 
maining objects, we centered the slit along the brightest 
nucleus identified in the optical/near-IR. 

All objects in the sample were observed in the Star- 
ing Mode Astronomical Observing Template (AOT) with 
the galaxies placed at two nod positions located approxi- 
mately one-third and two-thirds of the length of the slit. 
We were able to subtract out background sky contribu- 
tions in all modules for most targets; the CSOs observed 
in 2007 have no sky observations for the SH module. Ex- 
posure times were chosen to yield signal-to-noise ratios of 
S/N = 50 for the low-resolution modules and S/N > 20 
for the continuum in the high-resolution modules (Ta- 
ble©. 

The data were processed using the Spitzer Science Cen- 
ter S17.0 data pipeline. We used the basic calibrated 
data products, already corrected for flat-helding, stray 
light contributions, non-linear responsivity in the pix- 
els, and "drooping" (increases in detector pixel voltage 
during non-destructive readouts). The 2-D images were 
first medianed at each nod position to remove transient 
effects such as cosmic rays. For the SL and LL modules, 



we subtracted the sky contribution by differencing the 
images for each nod with the adjacent position in the 
same module. 

Following the initial cleaning, further rogue pixels were 
eliminated using IRSCLEAN.MASK0. We used rogue 
pixel masks for each IRS campaign and supplemented 
the standard masks with manual cleaning of each nod 
and module. The 1-D spectra were then extracted using 
the Spitzer IRS Custom Extractor (SPICE) v.2.0. We 
used the optimal extraction routine with the standard 
aperture to improve the S/N ratio. 

The low-resolution modules were stitched together to 
match their continuum levels by using a multiplicative 
scaling. The flux level of the LL1 module was kept fixed, 
since this module has the most overlap with the trans- 
mission filter of the 22 fim peakup. We then scaled the 
LL2 module to LL1, SL1 to the scaled LL2, and SL2 
to the scaled SL1. The average modular scaling factors 
were 0.98 for LL2 to LL1, 1.11 for SL1 to LL2, and 1.16 
for SL2 to SL1. We then calibrated the low- resolution 
spectra as a single unit by matching the convolved flux 
to the IRS 22 urn ded icated peakup using the method 
of lArmus et~aTI (|2007j ). For 4C +12.50, we scaled the 
IRS data to the IRAS 25 /im flux (10% uncertainty) 
with a scaling factor of 1.28. The required scaling in 

3 IRSCLEANJVIASK is available at 
http : //ssc . spitzer . caltech . edu/archanaly/ 
contributed/irscle WlRSCLEAN_MASK.html 
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the majority of cases was quite small, indicating that 
the sky subtraction and spectral extraction techniques 
are robust; the mean overall scaling factor for the IRS 
peakups was 1.04. Final spectra were produced from an 
error-weighted mean of the two nods. 

Noisy areas on both ends of the SH and LH orders 
were trimmed from the final spectra. These areas typi- 
cally encompass a range of 10-30 pixels on the edge of 
the orders and correspond to areas of decreased sensi- 
tivity on the detector. We trimmed only pixels with an 
overlapping wavelength range in adjacent orders so that 
the maximum amount of information is preserved. 

3.1. 1245+ 676 non-detection 

One of the CSOs, 1245+676, lies at a distance of 
495 Mpc (the third-most distant galaxy in the sample) 
and was predicted to have a flux density of ~ 10 mJy 
near 12 /zm, based on a linear interpolation between ex- 
isting optical/near- IR and radio photometry. This same 
method yielded flux estimates within a factor of a few for 
the other three CSOs in the sample without IRAS de- 
tections. The SUR peakup of 1245+676, however, only 
yielded flux densities of 0.7 and 1.0 mJy at 16 and 22 /im. 
The source of the flux estimation error likely comes from 
the 2MASS near-IR photometry, which may sample near 
the peak of an older population of red stars. If this is the 
case, then the true near-IR/optical continuum levels are 
significantly lower, which leads to a steeper optical-radio 
index and a lower estimate for the flux at 12 /zm. 

The order of magnitude deficit in the expected flux 
means that the integration times for the IRS were too 
short, and the extracted spectra have no identifiable fea- 
tures (at 15 /im, the LR modules have S/N< 3 and the 
SH module has S/N< 1). The peakup measurements in 
Tableware regarded as valid, but all features in the 1-D 
spectra are limits (including the spectral indices, which 
are consistent with no slope). Therefore, we do not in- 
clude data from this galaxy in subsequent analyses. 

4. RESULTS 

We show the photometrically-scaled low-resolution 
spectra in Figure [TJ with individual modules stitched to- 
gether. The SH spectra are shown in Figure [5] and the 
LH spectra in Figure [3] While individual orders within 
the high-resolution modules were typically well-aligned 
in flux, the differences in calibration between the SH 
and LH modules are clearly apparent when matching the 
spectra; this is due to a combination of different slit sizes 
for the SH and LH modules (by a factor of ~ 4) and the 
lack of separate sky subtraction for the SH modules in 
some objects. For this reason, continuum levels between 
the SH and LH modules do not always match; the LH 
module is typically the more accurate measurement due 
to its removal of sky background in all galaxies. 

There still exists significant noise in the reduced spec- 
tra, most likely due to hot pixels or other instrumen- 
tal artifacts from the IRS that were not corrected dur- 
ing the cleaning process. This occurs most often in the 
LH spectra - an example resembling a wide emission 
feature occurs between the 11 th and 12 th orders (see 
PKS 1413+135 near X rest = 27 /im in Figure I3J). We 
mark the locations of genuine features on the spectra 
themselves, with all detections we regard as real listed 



TABLE 3 

Peakup photometry and mid-IR spectral indices 



Object 


16 /tm PU 


22 PU 


"15-6 


"30-20 




[mJy] 


[mJy] 






4C +31.04 


4.8 


6.3 


0.6 + 0.3 


1.8 + 0.4 


4C +37.11 


5.5 


9.6 


1.1 + 0.6 


1.8 + 0.3 


1146+59 


18 


24 


0.3 + 0.2 


1.7 + 0.3 


1245+676 


0.74 


1.0 






4C +12.50 






2.3 + 0.2 


2.0 + 0.3 


OQ 208 


234 


396 


1.4 + 0.1 


0.9 + 0.1 


PKS 1413+135 


35 


47 


1.0 + 0.1 


1.8 + 0.4 


PKS 1718-649 


49 


57 


1.2 + 0.2 


0.0 + 0.1 


1946+70 


4.4 


6.0 


1.0 + 0.2 


0.4 + 0.4 



Note. — Errors in peakup (PU) fluxes arc at the 15% level. 
1245+676 has no spectral index measurements (sec S|3.H . 

in Tables SI [5j and |5] All other features are considered 
spurious. 

4.1. Continuum 

We measured the slope of the continuum in log space 
to create sets of rest-frame spectral indices for all objects 
(Table[3|). The continuum is characterized in two parts: a 
near-IR index between 6 and 15 /im (015-6) and a mid- 
IR index between 20 and 30 /im (030-20)- All indices 
are measured in regions of line-free continuum, verified 
by visual examination for each object. The slope of the 
continuum is closely related to the shape of the under- 
lying dust spectrum dominated by reprocessed thermal 
ra diation from dust . IRS spectra from the larger sample 
of lHao et al.l {2007) show a distinct change in spectral 
slopes according to their optical classifications: quasars 
have the shallowest 030-20 and 015-6, with AGN (Seyfert 
1 and 2) slopes growing steeper and ULIRGs having the 
steepest slopes in both regimes. 

The CSOs have a wide range of spectral indices; 
ai5_6 ranges from 0.3 to 2.3, and 0130-20 from 0.0 to 2.0. 
The steepest 015-6 and 030-20 ar e both in 4C +12.50, 
which has indices consistent with a single power law 
for the entire mid-IR spectrum. 1146+59 has the shal- 
lowest oi5-6 index and is flat between 5 and 8 /zm. 
PKS 1718—649 has the sh allowest 030 1-20 , also consis- 
tent with a flat spectrum. IWu et all ()2009l ) show that 
Seyfert 2 galaxies have steeper 030-20 (1-53 + 0.84) than 
in Seyfert 1 (0.85 ± 0.61). The two Sy2 galaxies in the 
sample (4C +37.11 and 4C +12.50) have 030-20 of 1.8 
and 2.0, while the only Syl CSO in the sample (OQ 208) 
has 030 - 20 = 0-9; fully consistent with this trend. The 
average indices for the sample are 015-6= 1.1 ± 0.6 and 
"30-20= 1-3 ± 0.8. 

Figure [4] shows overlays of the LR spectra for all nine 
CSOs, with fluxes normalized in the feature-free region 
near 15 /im, along with a composite template medianed 
over each channel. The continuum emission of the indi- 
vidual galaxies varies most strongly longward of 15 /zm; 
while most galaxies have 030-20 consistent with the tem- 
plate spectral index of 1.0, 4C +12.50 and 4C +37.11 are 
both significantly above the average. The steepness of 
oi5_6 for 4C +12.50 also strongly stands out, with the 
normalized spectrum falling below the template and all 
other galaxies shortward of 9 /im. 

4.2. PAH features 
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Fig. 1. — Full low-resolution (LR) IRS spectra of the CSOs, with detections of PAHs, silicates and water ice features marked for each 
object. 
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Fig. 2. — Short-high (SH) spectra of the CSOs. Detected fine-structure atomic (bottom) and H2 (top) emission lines are marked on each 
spectrum. PAH features are also marked along the bottom, although the broader profiles are more easily seen in the LR modules. 
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Fig. 3. — Long-high (LH) spectra observed with the Spitzer IRS. Symbols arc the same as Fig. [2] Features that are not labeled (such as 
the peak /im for PKS 1413+135) are considered spurious. 
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Fig. 4. — Low-resolution spectra (grey) of the CSO sample, 
shifted to rest frame and normalized in flux at a line-free area at 
15 /im. The thick spectrum is a composite template made from the 
error-weighted average of the individual scans, interpolated over a 
wavelength grid of AA = 0.08 ^m. 

Polycyclic aromatic hydrocarbons (PAHs) are among 
the dominant features in mid-IR spectra of dusty galaxie s 
(jHelou et al.ll2000t iDesai et al.ll2007t iSmith et "al1l2007h ; 
the emission is thought to be caused by UV excitation of 
vibrational modes in large molecules that are re-radiated 
in the infrared. We measured the strength of the pri- 
mary PAH features using two different app roaches: the 
first fits a spline to the lo cal continuum (|Spoon et al.l 
I2007t iZakamska et al.ll2008l) and then integrates the to- 
tal flux after a baseline has been subtracted. The sec- 
ond method uses PAHFIT (jSmith et al.ll2007[ L an IDL 
routine which performs a multi-component global fit on 
LR spectra based on a template for galaxies with a wide 
range of star formation. The advantage of PAHFIT is 
that it takes into account significant flux from the broad 
wings of dust features that blend with continuum in the 
spline-fit method; this is especially critical in measuring 
the equivalent widths (EWs) of PAH features. In order 
to compare the C SO hosts to other datasets (e.g., the 
"fork" diagram in iSpoon et al.l l2007l ). however, we fol- 
low the majority of the literature in using the spline-fit 
method. Results from both approaches are listed in Ta- 
ble n 

The strengths of the PAH features in the CSOs vary 
considerably, with no uniform template typical for the 
sample (Figure 13}. PKS 1413+135 shows no evidence 
for PAHs in any of the five main bands commonly seen 
in ULIRGs (6.2, 7.7, 8.6, 11.3, and 12.7 /im). This may 
be an indication of the hard UV/X-ray e nvironment that 
is known to dissociate PAHs (Voit 1992), which can arise 
from either extreme star formation or direct exposure to 
a central power source such as an AGN. It could also 
signify a lack of gas/dust near the ionizing source. No 
other tracers of star formation (such as Ly-a or fine- 
structure lines in the optical/IR bands) are identified in 
PKS 1413+135, suggesting that either the radiation of 
the AGN or large amounts of obscuration are responsible 
for the lack of PAH detections in this object. 

4C +12.50 and OQ 208 both show weak PAH emission 
above a strong continuum. The 11.3 and 12.7 /im fea- 
tures are strongest in OQ 208, while the PAH emission 



in 4C +12.50 is very dim for all detections. PAH fea- 
tures are seen in all five main bands for both objects, in 
addition to bands at 14.2 and 16.4 /im in OQ 208. We 
note that PKS 1413+135, OQ 208, and 4C +12.50 are all 
strong X-ray sources with L2-10 keV > 10 42 erg/s, fur- 
ther suggesting that dissociation of the molecules may be 
the primary reason for weak emission in these galaxies. 

1146+59, 4C +31.04, 1946+70, PKS 1718-649, and 
4C +37.11 all show moderate PAH emission strongest at 
11.3 and 12.7 /im. Two CSOs, 1946+70 and 4C +37.11, 
have no detectable 6.2 //m PAH emission, although the 
upper limit is consi stent with the 6.2/ 11.3 /im band ra- 
tio of 1.1 found by ISmith et al.l (|2007l ) for star-forming 
galaxies. Previous studies of Galactic PAH bands show 
that the 6.2 /jm luminosity can be suppressed with re- 
spect to the other bands for multiple reasons, includ- 
ing differential extinction ([Rigopoulou et a l. 2002) or hot 
dust increasing the continu um level and obs curing the 
wmgs of the PAH feature (jDesai et all 12007). A neu- 
tral ISM may also affect their relative strengths, as the 
6.2 uto. mode is stron ger when the PAHs are ionized 
(jVerstraete et al.lll996h . 

4.3. Silicate absorption 

Another prominent feature in the LR spectra is the 
pair of amorphous silicate emission/absorption com- 
plexes centered at 9.7 and 18 /im. The strengths of 
the silicate feature s are measured using the method of 
ISpoon et all (|2007l ): 

Sx = ^ (^A , (1) 

\ ± cont / 

where Fx is the flux extremum near 9.7 or 18 /zm and 
F CO nt is the continuum flux interpolated at the same 
wavelength from a spline fit. Results for the silicate 
strength fall into two categories: emission and weak ab- 
sorption (Table S|) . 

OQ 208 is the only object that shows silicate in emis- 
sion, with 5*9.7 = 0.3 and Sis = 0.2. The remainder of 
its mid-IR emission fo llows a relatively flat power law 
with an index of 1.9. lHao et al.l (|2007n found that sil- 
icate emission in galaxies occurred almost exclusively 
for quasars and Seyfert 1 galaxies, with an average 
59.7 = 0.20 for quasars and 5g.7 = —0.18 for Seyfert 1 
galaxies. The LR spectrum for OQ 208 closely resembles 
the average mi d-IR quasar temp late generated from ^20 
objects in the lHao et all (|2007| ) sample; this suggests a 
relatively small a mount of dust that is optically and ge- 
ometrically thin (|Levenson et al.l 120071 ) to produce the 
emission feature. A clumpy medium has been shown to 
reproduce such features due to illuminated clouds behind 
the central source (relative to the observer) filling in the 
absorption trough s ufficiently to create silicate emission 
(|Nenkova et al.ll2002t) . 

Six CSOs show weak absorption, all with silicate fea- 
tures ranging from 0.0 < 59.7 < —0.8. Narrow line emis- 
sion from H2 S(3) and/or [S IV] can be seen on top of the 
silicate absorption in several objects, suggesting that at 
least some part of the narrow-line gas is foreground to the 
absorbing dust. The average silicate strength is consis- 
tent with that se en both in Seyfe rt galaxies and ULIRGs 
(|Hao et al.l 120071) . Models from iLevenson et al.l ([20071 ) 
suggest that this type of emission may arise from dust 
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TABLE 4 

PAH AND SILICATE FEATURES IN LOW-RESOLUTION SPECTRA 





PAHFIT luminosity 


PAHFIT EW 


Spline-fit luminosity 


Spl 


nc-fit EW 


Silicate 


Object 


6.2 


11.3 


6.2 


6.2 ice 


11.3 


6.2 


11.3 


6.2 


6.2 ice 


11.3 


Sg.7 


Sl8 




[log L/L@] 


[log L/L Q ] 


[fJJJ\] 


[fJHx] 


H 


[log L/Lq] 


[log L/L@] 


H 


H 


H 






4C +31.04 


8.02 


8.19 


0.18 


0.16 


0.62 


7.88 


7.81 


0.18 


0.12 


0.32 


-0.8 


-0.4 


4C +37.11 




7.56 






0.28 


< 7.21 


7.25 


< 0.08 




0.17 


-0.5 


-0.0 


1146+59 


6.61 


7.01 


0.05 




0.37 


6.00 


6.71 


0.01 




0.23 


-0.7 


-0.2 


4C +12.50 


9.37 


9.43 


0.07 




0.05 


8.20 


8.73 


0.01 




0.01 


-0.5 


-0.0 


OQ 208 


9.24 


9.36 


0.08 




0.09 


8.75 


8.69 


0.03 




0.03 


+0.3 


+0.2 


PKS 1413+135 




< 8.82 






< 0.02 


< 9.05 


< 8.76 


< 0.01 




< 0.01 


-1.0 


-0.3 


PKS 1718-649 


7.11 


7.49 


0.08 




0.30 


6.67 


7.08 


0.04 




0.13 


-0.4 


+0.3 


1946+70 




8.04 






0.22 


< 7.48 


7.58 


< 0.04 




0.09 


-0.6 


-0.3 



Note. — The 6.2 )im PAH "ice" data use a continuum that has been corrected for water ice absorption (where present) at 6 fim. — indicates 
that PAHFIT fit no significant flux for a particular dust component. Silicate strength is defined in Equation [Tl 



clouds that are optically but not geometrically thick (ie, a 
clumpy ISM, characteristic of optically identified AGN). 
None of these galaxies show strong absorption from the 
secondary silicate feature at 18 /xm. PKS 1413+135 
shows the deepest silicate absorption at S9.7 — —1.0. No 
objects show the deep absorption features (Sg.7 < —1.5) 
seen in Seyfert 2 galaxies and ULIRGs that can only be 
generated with a smooth, spherical shell of dust. 

The 9.7 /im silicate depth is also translated into an ap- 
parent visual ext i nction using the Galactic calibration of 
iRoche fe Aitkenl (|T984); the measured extinction spans 
A v ~ 10-30 mags (Tabled]). 

4.4. Other absorption features: water ice, hydrocarbons, 
and gas-phase molecules 

Water ice absorption at 6 /xm has been observed in the 
mid-IR in ULIRGs but is typically absent in Seyfert and 
starburst galaxies ([Spoon et al.l l2000l l2002t ). We detect 
water ice in 4C +31.04 at an optical depth of r = 0.57. 
The feature is shallow and bracketed by possible PAH 
emission near 6.2 /xm, which makes accurate measure- 
ments of the local continuum difficult. There is no indi- 
cation of water ice in any of the other objects, nor in the 
stacked LR spectra of the CSOs (Figure H]); this argues 
against the existence of any weak ice features not seen in 
individual objects. 

Bending m odes in hydrogena t ed amorphou s car- 
bon (HAC) (iSpoon et all JpOOH iDartois et all [20071: 
IDartois fc Mufioz-Cardl2007f) cause two blended absorp- 
tion features in the mid-IR centered at 6.85 and 7.25 /xm. 
We detect both features at a weak level in 1146+59; 
we use a spline-fit to characterize the local continuum 
and measure their optical depths and integrated fluxes. 
The 6.85 /xm HAC has r = 0.12 and / = -0.7 x 
10~ 21 W cm" 2 , with the 7.25 /xm HAC at r = 0.10 and 
/ = -0.6 x 10~ 21 W cm" 2 . The CSO detection rate for 
HACs (1/8) is consistent with that seen in larger samples 
of IR-bright galaxies (Willett et al. 2010), especially in 
those where strong emission from the 6.2 /xm PAH fea- 
ture may obscure the HAC absorption. 

Gas-phase molecules (including C2H2, HCN, and CO2) 
have also been observed in absorption in galaxies in 
the mid-IR, typically in deeply obscured ULIRG nu- 
clei (|Lahuis et al.l 12007'). No detections of any of these 
molecules were made for the sample; given the character- 
istic column density needed to show absorption, however, 
the S/N is likely too low in all galaxies to be detected at 



a typical column density of N > 10 16 cm 2 . 

4.5. Atomic lines 

The high-resolution data is primarily used for the iden- 
tification of atomic and molecular emission lines with 
widths too narrow to be separated in the LR modules. 
Integrated luminosities (assuming isotropy) are given in 
Table O upper limits are computed by assuming a Gaus- 
sian with a peak flux density 3cr above the noise and a 
velocity width derived from detections in the sample of 
ULIRGs from Willett et al. (2010). Line fluxes were 
measured using the standard packages in the Spectro- 
scopic Modeling Analysis a nd Reduction Tool (SMART) 
v6.2.4 ([Higdon et all 12004) . We found a simple Gaus- 
sian to be a good fit for virtually all high-resolution lines, 
with the central wavelength left as a free parameter in 
the fit; for cases in which lines are blended (such as the 
[O IV] /[Fe II] complex), we used a two-component Gaus- 
sian fit centered at their expected wavelengths from op- 
tical redshifts. 

The fine structure [Ne II] 12.8/xm and 
[Ne III] 15.6/xm transitions are the most commonly 
detected in the sample, with both observed in seven 
objects. We also detected [Ar III] 9.0/xm, [S III] 18/xm, 
S III] 33/xm, [S IV] 10.5/xm, [Si II] 35/xm, [Fe II] 23/xm, 
Fe II] 26/xm, and [O IV] 25.9/xm. The strongest lines 
seen in the HR modules, [Ne II] and [Ne III], were 
resolved in 6/7 objects with neon detections. The 
velocity width for these transitions spans 250-600 km/s, 
with the broadest example the [Ne III] transition in 
4C +12.50 (see ^5.51 for a discussion of linewidths for 
this galaxy). The high- ionization [Nc V] line, commonly 
observed in AGN with accreting gas, was not detected in 
any object in either the 14.3 or the 24.3 /xm transition. 

As a check for our line measurements, we compared 
lines detected in the arc hival galaxy 4C + 12.50 with 
HR results published by IFarrah et alj (|2007f) . For the 
most part, the two sets of measured fluxes are in good 
agreement; however, we cannot confirm detections of the 
[Ne V] 14/xm or H2 S(0) transitions, both o f which are 
considered uncertain by IFarrah et alj |2007). The au- 
thors also measure a flux for the [Ne III] transition that 
is higher than our measurement by ~ 30%; this may be a 
result of different fitting techniques, since the line has a 
non-Gaussian component in a broad blue wing, p ossibly 
from a photoionized outflow ([Spoon fc Holt 200J). 
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TABLE 5 

Atomic line luminosities from HR spectra 



Object 
A res t \pm] 


[Ar m] 
8.991 


[S iv] 
10.511 


[Ne n] 
12.814 


[Ne v] 
14.322 


[Ne in] 
15.555 


[Fe n] 
17.936 


[S in] 
18.713 


[Ne v] 
24.318 


[O iv] 
25.890 


[Fe n] 
25.988 


[S m] 
33.481 


[Si n] 
34.815 


4C +31.04 




< 7.03 


7.34 


< 6.89 


7.32 


< 6.53 


6.71 


< 6.72 


7.04 


6.61 


< 7.39 




4C +37.11 




< 6.71 


7.81 


< 6.63 


7.43 


< 6.38 


6.97 


< 6.59 


6.43 


6.93 


7.30 


7.93 


1146+59 




< 5.72 


6.19 


< 5.51 


6.12 


5.40 


5.64 


< 5.24 


5.95 


5.18 


6.11 


6.34 


4C +12.50 


7.76 


8.19 


8.66 


< 8.01 


8.62 


< 8.02 


8.03 


< 8.33 


8.55 


< 8.55 






OQ 208 




< 7.79 


8.22 


< 7.66 


8.10 


< 7.39 


7.68 


< 7.54 


7.47 


7.45 


< 8.33 




PKS 1413+135 


< 8.34 


< 8.58 


< 8.44 


< 8.27 


< 8.76 


< 8.23 


< 8.35 


< 8.20 


< 8.49 


< 8.54 






PKS 1718-649 




< 5.97 


6.80 


< 5.66 


6.58 


< 5.47 


6.39 


< 5.55 


5.75 


5.87 


6.35 


6.95 


1946+70 


< 6.73 


< 7.22 


7.58 


< 6.74 


7.35 


< 6.94 


< 7.32 


< 7.41 


< 7.56 


< 7.62 







Note. — All line luminosities are in log Lq- : — indicates that the observed wavelength was not visible in the IRS range. Errors in 
the flux measurements are on the order of 10%. 



4.6. Molecular hydrogen 

Lines from the rotational series of H2 are seen in 7/8 
CSOs. The exception, PKS 1413+135, showed no iden- 
tified atomic or molecular emission lines at all in the 
mid-IR. The strongest line in all objects is the ortho 
H2 S(l) 17.0/xm transition, appearing in every object in 
which molecular hydrogen was detected. The next ortho 
transition, S(3) 9.67/xm, was detected in five objects; we 
can only place upper limits on the S(5) 6.91/^m transi- 
tion due to likely blending from [Ar II] 6.99/xm in the 
low-resolution spectra. For the para H2 transitions, the 
ground state S(0) 28.2/xm was detected in only one ob- 
ject, although this line can be obscured by rising con- 
tinuum levels near 28 /im. The next excited state, 
S(2) 12.3/xm, is detected in six objects, with OQ 208 the 
only object in which the S(l) transition was observed 
without S(2). Fluxes for all lines are given in Table [6j 
upper limits are computed in the same manner as those 
for atomic emission lines. 

4C +37.11 shows a highly unusual H2 spectrum, with 
all lines from S(l) through S(7) detected in the LR data. 
This is the only object that shows clear detections of the 
S(4) 8.03/xm, S(6) 6.11/zm, and S(7) 5.51/im lines, and 
the H2 luminosities are much larger compared to the 
strength of its fine-structure lines (such as neon) than 
in any o ther object. I n a s ample of 77 ULIRGs ana- 
lyzed bv lHigdon et al.l (|2004D . the S(4)-S(7) lines were 
seen only in a single nearby target (NGC 6240). Sim- 
ilar H2 strengt hs are also obse rved in the FR II radio 
galaxy 3C 326 (O de et al~ll2007l) and the compact group 
Stephan's Quintet (|Appleton et al.l [2006) . Since other 
tracers such as PAH strength and [Ne II] emission are in- 
consistent with extreme amounts of star formation in this 
galaxy, the strong H2 lines in these objects are attributed 
to shock heating of the gas by a tidal accre tion flow aris- 
ing fro m interactions with other galaxies. iManess et al.1 
(2004) show that the radio structure in 4C +37.11 may 
come from interactions driven by a recent merger, with 
a dense circumnuclear medium giving rise to shocks in 
the central component that would be consistent with the 
H2 lines seen here. 

5. INFRARED DIAGNOSTICS 

5.1. Dust and gas temperatures 

We used mid-IR phot ometry from the IRS peakups and 
published IRAS fluxes dSanders et al ll2003D to fit a mod- 
ified blackbody curve (|Yun fe Carilhll2002l ) which fits for 
a dust temperature (Tdust) and an angular size (Table [7]). 



Since the CSOs observed in 2009 were undetected by the 
IRAS satellite, fits for these objects are based only on the 
two IRS peakup point s and are too unders ampled for a 
measurement of Tdust ■ lYun fe Carillil (|2002T ) estimate the 
scatter on Tdust in their sample to be ~ 10K, while the 
formal errors on the curve-fitting give errors on the order 
of a few K. 

Out of the four CSOs with IRAS fluxes, 4C +12.50, 
4C +31.04, and OQ 208 are moderately well-fit by a 
single-temperature blackbody, with Tdust ranging from 
65-110 K. The data for PKS 1413+135 has signifi- 
cant flux excesses to the fitted curve at higher ener- 
gies; possible causes include contributions from hotter 
dust (500-2000 K), significant line emission in the 12 
and 16 /im bands, or non-thermal emission. This last 
option is supported by identification of PKS 1413+135 
as a BL Lac object, with beaming in the core con- 
tribut ing to continuum- dominated NIR and MIR spec- 
tra (jPerlman et al.lll996| ). The blackbody fit for 1146+59 
shows flux excess above the fitted curve at lower energies, 
which may indicate large amounts of dust at cooler tem- 
peratures (< 50 K). As a sample, the relatively poor fits 
to a modified blackbody most likely indicate that these 
galaxies ar e complex systems without a single character- 
istic Tdust- iDeo et al. (2009|) show that Seyfert galaxies 
have at least three thermal components at T ~ 1000, 
200, and 60 K. 

For the galaxies with multiple transitions of molecu- 
lar hydrogen, we also fit excitation temperatures (T ex ) 
to a warm H2 compo nent following the meth ods of 
iRigopoulou et all (|2002t l and lHigdon et al.1 (|2006f ). Tem- 
peratures are calculated by assuming a Boltzmann distri- 
bution for the series of lines, the ratios of which depend 
only on T ex ; galaxies in which lines higher than S(3) are 
detected are typically better fit with a two-temperature 
model. The CSOs span a broad range of tempera- 
tures, with T ex ranging from 202 K in PKS 1718-649 to 
537 K in 1146+59 (Figured]). The line luminosities can 
also be used to measure masses of the H2 components 
(Table [7]) ; the mass of the warm gas varies over more 
than two orders of magnitude for the CSOs, with the 
smallest contribution coming from the gas in 1146+59 
{Mh 2 = 9 x 10 4 M ) to the much larger component in 
the galaxy 4C +12.50 (M a = 4 x 10 7 M Q ). 

4C +37.11 is the only galaxy for which a hot gas tem- 
perature can be measured, at T = 854 K; the derived 
hot gas mass is ~ 10% of the warm H2 gas. Measure- 
ments of H 2 from lHigdon et all (|2006l ) and lRoussel et al.l 
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TABLE 6 

Line luminosities of molecular H2 gas 



Object 


H 2 


S(7) 


H 2 


S(6) 


H 2 S(5) 


H 2 


S(4) 


H 2 


S(3) 


H 2 


S(2) 


H 2 S(l) 


H 2 


S(0) 


A reS i \pem] 


5 


51 


6 


11 


6.91 


8 


03 


9 


67 


12 


.28 


17.04 


28 


.22 


4C +31.04 


< 


7.34 


< 


7.25 


7.53: 


< 


6.47 




7.37 




6.87 


7.46 


< 


6.77 


4C +37.11 




7.07 




6.66 


7.70: 




7.09 




7.59 




7.18 


7.38 


< 


6.98 


1146+59 


< 


6.09 


< 


5.91 


6.08: 


< 


5.23 


< 


5.29 




5.74 


5.75 


< 


5.49 


4C +12.50 


< 


8.53 


< 


8.38 


< 8.84 


< 


7.88 




8.21 




8.04 


8.43 


< 


8.55 


OQ 208 


< 


8.07 


< 


7.93 


8.27: 


< 


7.49 




7.31 


< 


7.50 


7.86 


< 


7.51 


PKS 1413+135 


< 


8.80 


< 


8.62 


< 8.79 


< 


8.24 


< 


8.38 


< 


8.41 


< 8.22 


< 


8.79 


PKS 1718-649 


< 


6.41 


< 


6.24 


6.69: 


< 


5.57 


< 


5.83 




6.49 


6.87 




6.26 


1946+70 


< 


7.26 


< 


7.08 


< 7.53 


< 


6.58 




7.44 




7.16 


7.67 


< 


7.57 



Note. — Luminosities arc given in log Lq. Measured H2 S(5) lines (indicated with :) are upper limits 
due to possible blending with [Ar 11] emission; sec jj4,6| for details. The H2 S(4-7) lines were measured in the 
LR modules; all other lines were measured in the high-resolution modules. Errors in the flux measurements 
are on the order of 10%. 



TABLE 7 

Derived mid-IR. properties of the CSO sample 



Object 


SFR]V e 


SFRp AH 


A v 


rpwarvn 




rphot 




^ dust 


log M°™ 


log M b B ^ e 




[M /yr] 


[M e 


M 


[mag] 


[K] 


[1O 7 M ] 


[K] 


[1O 7 M ] 


[K] 


[Mq] 


[Mq] 


4C +31.04 


7.8 + 1.1 




6.4 


15 + 6 


338 ± 100 


0.47 + 0.13 






65 


< 8.16 


8.78 


4C +37.11 


17 + 1 


[0.8 - 


-i.fi] 


9 + 3 


354 ± 90 


0.95 + 0.05 


854 ± 240 


0.09 + 0.02 




< 7.96 


8.74 


1146+59 


0.5 + 0.1 


0.3 


13 + 1 


537 ± 50 


0.01 + 0.002 






150 


9.23 + 0.52 


8.40 


4C +12.50 


159 ± 14 




31.5 


8 + 0.6 


326 ± 97 


4.26 + 0.46 






85 




8.81 


OQ 208 


54 + 3 




47.8 




281 ± 101 


1.13 + 0.13 






110 


< 7.79 


8.81 


PKS 1413+135 


< 150 




< 80 


13 + 0.4 










90 




8.22 


PKS 1718-649 


1.8 + 0.1 




0.8 


7 + 3 


202 ± 112 


0.12 + 0.01 








8.62 ± 0.45 


8.48 


1946+70 


11.1 + 0.9 


[1.7- 


-3.1] 


11 + 3 


320 ± 83 


0.74 + 0.05 










8.54 



Note. — Multiple values for SFRpah indicate that only the 11.3 /^m PAH feature was detected, giving a range in possible star formation rates. 



BH masses are calculated both from the [OIV] lincwidth and the L 



bulge 



■ M B h relation of Bcntz ct al. (2009). Errors i 



M 



bulge 



arc ~ 0.1 dcx. 



(2007) combined with CO measurements suggest that the 
warm/hot components probed by the mid-IR lines are 
typically a small fraction of the total molecular gas mass 
(1 — 30%), the majority of which exists at much lower 
temperatures. 

The excitation temperatures measured for the warm 
H2 gas are typically a factor of several higher than the 
dust t emperatures measur ed by fitting the IR photom- 
etry. iHigdon et al.l ()2006l ) found no tendency for the 
H2 temperatures to depend on dust temperatures, based 
on IRAS fluxes. Results from the SINGS sample suggest 
that much of the H2 emission comes from either photodis- 
sociation regions (PDRs) associated with Hll regions in 
star-forming galaxies, or from nuclear regions close to 
low-level AGN with an additional H2 component from 
shock heating (jRoussel et al.1 120071 ). In both cases, the 
gas and dust are not expected to be co-spatial, and there- 
fore temperature differences in a multi-phase medium arc 
not unexpected. 

5.2. Silicates 

ISirockv et al.l (|2008f ) recently showed that the depths of 
the 9.7 and 18 /im silicate features depend sensitively on 
the geometric distribution of the dust. They identify two 
distinct classes of ULIRGs, seen both in IRS data and 
in numerical models that distinguish between dust that 
is clumpy (AGN with broad optical lines) and smooth 
(geometrically and optically thick, with lower amounts 
of ionizing radiation). Plotting the relative strengths of 
the dust features against each other yields a "feature- 



feature" diagram (Figure[6]), characterized by two regions 
representing the clumpy and smooth families in the op- 
tically thick limit. All Type 1 AGN in the sample of 
ISirockv et all (pOOrl are found only in the clumpy region; 
LINERs and Hll regions generally show deeper absorp- 
tion at 9.7 /Ltm, lying along the smooth track. Type 2 
AGN may belong to either family. 

Seven of the CSOs lie along the track for clumpy dust. 
Six of these galaxies have silicate ratios with A = 3 — 5 
(the average number of dust clouds along a line of sight). 
The other galaxy lying in the clumpy region, OQ 208, is 
the only object in the sample showing the 9.7 ^m feature 
in emission, indicating that the dust is optically thin. Its 
location in Figure [6] indicates that it has only No ~ 1 — 2. 
In the standard AGN model, this could be explained by 
viewing a central power source and associated dust torus 
close to face-on, with a high inclination angle reducing 
the possibilities for absorption along the line of sight to 
the observer. This agrees with the optical classification 
of OQ 208 as a broad-line radio galaxy. 

T he silicate rati o for P KS 1718—649 is consistent with 
the ISirockv et al.l (|2008l ) models for smooth, geometri- 
cally thick spherical shells of dust. The galaxy has an 
opt ical depth of ' ry ~ 100; there are few ULIRGs from 
the iSpoon et all (|2007[ ) and llmanishi et al.1 ([2007) sam- 
ples nearby in the feature-feature diagram. The host 
galaxy has an "incomplete" Hi ring and trailing spiral 
arms, both of which argue against a smooth distribution. 
If so, this may indicate that the absorbing dust exists on 
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Fig. 5. — H2 excitation diagrams for the CSO sample. Dotted lines show the T ex fit to the warm H2 gas, while dashed lines show the fit 
to a hot component (where detected). S(5) has error bars when the blended H2/[Ar II] feature is detected, but has only an arrow when it 
is a true upper limit. 

smaller scales (less than 1 kpc) close to the center of the 
galaxy. 

We tested the robustness of the feature-feature dia- 
gram by running numerical models on the CSOs to quan- 
tify their physical parameters. We used the CLUMPY 
code to model a toroidal distri bution of point-like d ust 
clouds around a central source (jNenkova et al.ll2008allbT ) . 
Individual clouds have optical depth ry and an expo- 
nential distribution such that the number of clouds per 
radial ray (No) is described by a Gaussian with angu- 
lar width cr in the polar direction. The thickness of the 
torus, Y , is scaled to the inner radius set by the dust 
sublimation temperature (Y = Router/ Rinner)- The fi- 
nal input parameters are the inclination of the disk i and 
the power law index q for the radial cloud distribution. 
The spectra are dominated by reprocessed dust emission 
independent of the heating source, but the AGN emis- 
sion is visible for models with low optical depths and 
high inclination angles. 

For the CSOs, we sampled a coarsely-defined grid of 
40,000 models across all six parameters. To assess the 
goodness-of-fit, we rebinned the IRS LR spectrum to the 
same wavelengths as the CLUMPY m odel and minimize 
the error E from lNikutta et al.l (|2009| ). This fitting also 
includes a multiplicative scaling of the spectrum account- 
ing for the total bolometric flux of the AGN. 

Since CLUMPY only models the emission from black- 
body dust, the fit is poor if the galaxy has strong emission 
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Fig. 6. — Feature-feature diagram displaying the relative 
strengths of the 9.7 and 18 fivci silicate features. The lines represent 
modeling results o f different dust geo metries for cool, oxygen-rich 
silicates (OHMc) HSirockv et alj|200gfl . Dotted, dashed, and solid 
lines represent the thickness of the dust shell in the spherical mod- 
els (Y = Router I Rinner) , and q is the power law index for the den- 
sity distribution of the shell (p[r] oc r~ q ). The green tracks show 
different numbers of dust clouds located along the line of sight for 
the clumpy model (No = 1, 3, and 5 from upper right to lower 
left). The black star in the upper right corner is the starting point 
for all optically thin dust in the models (Sg.7 = 1.26, Sis = 0.67). 
The CSOs are shown in black circles with error bars (smaller than 
the symbol in several cases). 
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Fig. 7. — IRS low-resolution spectra of three CSOs (top to bot- 
tom: OQ 208, PKS 1413+135, and 1946+70) with the best-fit 
models from the CLUMPY code overlaid in red. The minimized 
errors for the spectral fits are E m i n = 1.52, 0.50, and 0.20, re- 
spectively. Parameters of the CLUMPY fits are discussed in SI5.2I 



features from atomic, molecular, or PAH features. In the 
CSO sample, the best fit models have E m i n between 0.20 
and 1.52 - we show examples of three fits to the CSOs 
in Figure The best fit is for 1946+70, with Y = 100, 
A = 10, q = 1.0, and t v = 10. The model reproduces 
absorption in both the 10 and 18 fim features, but under- 
predicts emission shortward of 6 /im. This is a general 
trend for most CSOs in our sample, possibly due to hot 
dust (T ~ 1000 K) or PAH emission. Neither a nor i are 
well-constrained for any set of models. 

Two other CSOs show interesting fits from the 
CLUMPY models. OQ 208 exhibits the worst formal fit 
in the sample, with E m i n = 1.52. This is almost entirely 
due, however, to the underprediction of flux shortward of 
9 fxm by ~ 0.1 dex. The model fits the spectrum well be- 
tween 9-28 /im, tracing both silicate features in emission. 
The best-fit models for this galaxy favor Ao between 5 
and 20, q = 2.0, t v between 200 and 300, and Y be- 
tween 30 and 100. These parameters agree with general 
predictions of iNenkova et al . (2008b) and its placement 
in Figure |5] with the exception of an unusually high ry > 
rarely seen in Seyfert galaxies. 

PKS 1413+135 is another galaxy of particular interest 
due to its unusual spectrum - the lack of any atomic or 
molecular features make it an ideal candidate for radia- 
tive transfer modeling. The best-fit CLUMPY model for 
PKS 1413+135 has Y = 10, A = 10 - 20, q = 0.0 - 1.0, 
and Ty = 10 — 30. The model fits both dust absorption 
features well, but severely underpredicts the flux on ei- 
ther side of the silicate features (more than 0.5 dex at 
5 and 30 /im). Both q and Ty agree with its placement 
in Figure [51 but the torus thickness Y is unusually thin. 
The excess at the blue end could be due to additional 
reservoirs of hot dust not modeled by CLUMPY; how- 
ever, the blue excess would likely require large amounts 
of very cold dust not supported by the amount of star 
formation seen. The best-fit model of PKS 1413+135 
with a smooth spherical geometry yields Y — 100 — 500, 
q = 1.5, and Ty = 60 — 90; this fits the blue excess well, 
but fails to account for the 18 /jm absorption and the 
red excess. The optical depth values are similar to values 



-5 
-4 

_c 

o> -3 
c 

CD 

cn _ 2 

CD 
O 

HE -1 



ty) 







i — i — i — 

CSOs • 
AGN 

Obscured AGN X 
Normal ■ 
Storbursts ♦ 
ULIRG/HyLIRGs * 



* « 



0.001 0.010 0.100 1.000 

6.2 urn PAH EW |>m] 

Fig. 8. — "Fork" diagram displaying the 6.2 /im PAH equivalent 
width vs. the silicate strength at 9.7 /im. CSO galaxies are s hown 
in red; additional Spitzer data are from [Spoon ct al. ( 2003) an d 
Willett et al. (2010). Limits are given on the PAH EW for CSOs 
with no measured flux above the noise. The black arrow gives the 
mean effect on the PAH EW if fluxes from PAHFIT are used in 
place of the spline-fit technique. 

obtained from both the X-ray spectrum (P erlman et al., 
in pr e p) and the Hi 2 1 -cm a bsorption line (Caril lTet al.l 
1992). iPerlman et al. ( 2002) propose a geometry with a 
single CMC lying along the line of sight, which agrees 
with the low Y seen in our model. 

The individual CLUMPY models broadly support 
the conclusions from the feature-feature diagram. A 
clumpy dust geometry is a reasonable fit for most mod- 
els, although there are clear contributions from dust 
at cold temperatures that are not properly modeled. 
PKS 1718—649, the only galaxy predicted to have a 
smooth dust shell, does not give a good fit to either sili- 
cate absorption feature with the CLUMPY geometry. 

5.3. Quantifying the AGN fraction 

A powerful mid-IR diagnostic for separating popula- 
tions of starburst g alaxies from AGNs is the "fork dia- 
gram" developed bv lSpoon et al.l (|2007D . which plots the 
6.2 /im PAH equivalent width (EW) against the depth 
of the 9.7 /jm silicate absorption. Two distinct families 
are seen on this plot, with a horizontal branch of low 
silicate absorption/emission and variable PAH content 
largely populated by AGN, and a diagonal branch with 
deeper dust absorption and stronger PAH emission pop- 
ulated mostly by starburst galaxies. We overlay the data 
from the CSOs on the fork diagram containing a larger 
sample of galaxies in Figure [H 

4C +31.04 lies between the two branches, with a PAH 
EW that does not clearly classify it as either a starburst 
galaxy or an AGN. The remaining CSOs in the sample lie 
along the horizontal branch in the locus of AGN galaxies, 
while four have upper limits on the PAH EW that may 
place them even further from the starburst branch. 

The separation of galaxies into two distinct branches 
is related to the geometric distribution of the dust, as 
shown in Figure [BJ The horizontal branch in Figure [5] is 
defined by the low amounts of 9.7 fim absorption found in 
clumpy models, with the varying PAH EW a function of 
either the ionizing flux from an AGN (which dissociates 
PAHs) or the metallicity of the galaxy. The diagonal 
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Fig. 9. — The PAH 6.2 (im equivalent width vs. the flux ratio of 
the [O IV] 25.9 (im and [Ne II] 12.8 fim transitions. We compare 
the CSOs to sample s of AGN, ULIRGs, and starburst galaxies from 
lArmus et all (200% . The mean of the AGN (red) and starburst 
galaxies (green) are assumed to set the 100% level of AGN and star 
formation contributions, respectively - the bars show the fractional 
contribution assuming simple linear mixing. 

branch contains the majority of the galaxies with deep 
silicate absorption and requires a geometrically smooth 
distribution; while this can harbor an AGN at the highest 
levels of obscuration, none of the CSOs are located in this 
region. This suggests that none of these CSOs have large 
amounts of star formation often seen in ULIRGs. 

A second method of assessing the AGN contribu- 
tion uses the 6.2 /im PAH EW plotted against the 
[O IV] /[Ne II] ratio (Figure EJ} . Since [O IV] is known 
to arise from both AGN and star formation, comparing 
the two addresses the source of these moder ate- ionization 
lines. For comparison, we use data from lArmus et al.l 
(2007|) where we assume that the [O IV] emission from a 
pure AGN represents 100% of the contribution to Lju; 
similarly, we use a sample of pure starburst galaxies to 
approximate the contribution from star formation. A 
simple linear mixing is assumed for tracing the distribu- 
tion between the extremes. 

Two CSOs (4C +12.50 and 1146+59) fall in the AGN 
family, with < 2% of the expected PAH contribution 
from star formation. PKS 1413+135 has no [O IV] or 
[Ne II] lines, but a PAH EW upper limit of < 0.1 /mi. 
The remainder show a larger star formation contribution 
to the [O IV] emission, with the highest starburst fraction 
at ~ 25%. Since no CSO shows either AGN or star 
formation contributions above 30% from this diagnostic, 
this indicates that CSOs form an intermediate class with 
substantial contributions from both energy sources. 

Finally, we compare the star formation contribution to 
the amount of hot dust present in the galaxies. AGN can 
heat silicate dust to temperatures in excess of 1000 K, in- 
creasi ng the continuum emission at shorter wave l ength s 
(e.g., iPier fc Krolikl [1991 lAlonso-Herrero et al.l |200TI ). 
Figure [TU] compares the PAH strength to the flux ratio at 
15 and 5.5 /im, which will be lower as the amount of hot 
dust increases. The three vertices represent "pure" ex- 
amples of an AGN, an Hll region, and a photo-dissociated 
region (PDR), along with data for AGN, starburst galax- 
ies, ULIRGs, and CSOs. PKS 1413+135 is the only CSO 
that falls strongly in the AGN region - the remainder of 



Fig. 10. — Peak PAH flux at 6.2 /mi (normalized to 
the 5.5 /im continuum) vs. the flux ratios at 15 and 
5.5 /im (Laurent et al. 2000). The PAH emission traces star for- 
mation, while the f(15)/f(5.5) ratio is inversely proportional to the 
amount of hot dust in the galaxy. Symbols are the same as in 
Fig. [9] The dotted line vertices represent "pure" examples of an 
AGN (3C 273), an Hll region (M17), and a photo-dissociation re- 
gion (PDR; NGC 7023); dashed lines indicate the fractional AGN 
contribution. 

the galaxies are scattered with ULIRGs indicating contri- 
bution from both star formation and AGN heating. The 
scatter between the two families is consistent with the 
trends seen from both the silicate depths and [O IV] ra- 
tios (Figs. |S]and[H]). There is no obvious correlation of 
the AGN fraction with silicate depth, kinematic ages of 
the jets, or optical type. 

The mix of AGN and star formation activity in all three 
diagnostic plots supports the picture that the active radio 
phase may be triggered following a merger. T he 10 8 yr 
delay seen in CSOs by iPerlman et al.l (|2001bD provides 
enough time for a burst of star formation triggered by the 
initial pass; afterwards, as gas is driven to the center of 
the merging galaxies, the AGN phase begins to dominate 
the IR emission. The CSO outburst, however, appears 
to take place well before the AGN overwhelms star for- 
mation as a photon source in the IR. The scatter in AGN 
fraction in the CSOs may have several sources, includ- 
ing differences in viewing geometry, the total amount of 
star-forming gas, and the time since the beginning of the 
active phase. 

5.4. Star formation rates 

We quant ify star formation rates (SFR) using the re- 
lation from IHo fc Ketol 1)20071 ) , where the intensities of 
mid-IR neon emission show a tight linear correlation with 
both Ljr and the ionizing Lyman-a continuum from 
young stars. The combined luminosity from both the 
[Ne II] and [Ne III] lines are used to determine the SFR: 

SFR [M /yr] = 4.73 x 10~ 41 [L NeII + L NeIII ] , (2) 

where th e neon luminositie s are measured in erg s _1 . 
Following LHo fc Ketol (|2007t ). we assume f ion = 0.6 (the 
fraction of ionizing photons trapped by the gas) and 
that the fractions of singly- and doubly-ionized neon 
are / + = 0.75 and / ++ = 0.1. We compare this cal- 
ibration for ULIRG data to the LiR-L^e fit for the 



MID-IR SPECTROSCOPY OF COMPACT SYMMETRIC OBJECTS 



15 





9.5 : 




9.0 r 


1 


8.5 r 


\ 
s 


8.0 \ 


_l 


7.5 r 


cn 




_0 


7.0 r 




6.5 \ 




6 ' 




9.0 9.5 10.010.511.011.512.012.5 
log [L IR /L Q ] 



Fig. 11. — Luminosity of [Ne II] + [Ne III] lines as a function of 
Lin f° r the four CSOs with IRAS detections. Solid lines show the 
best linear fit for the CSOs; the dotted line shows the fit for the 
much larger sample of Ho & Kcto (2007). An upper limit is given 
for PKS 1413+135, which shows no detection of either neon line. 
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Fig. 12. — Combined luminosity of the 6.2 and 11.3 /im PAH 
features as a function of the [Ne II] + [Ne III] luminosities. The 
dashed line s hows the linea r fit for the PAH-neon relation calibrated 
for ULIRGs (Farrah ct al. 2007); the solid line shows the linear fit 
to the CSOs with neon detections. 

CSOs (Figure [TTj) . IRAS fluxes are available for four 
CSOs in the sample, w hich we use to compute Lib. 
(|Sanders fc Mir abel 19961) : the fit is reasonably good for 
all objects, even with the CSOs varying in Lm by more 
t han three orders of m agnitude. 

iFarrah et all (|2007h use a second method of determin- 
ing star formation rates that makes use of the observed 
linear correlation in ULIRGs between the luminosities of 
[Ne II] + [Ne III] and the 6.2 and 11.3 fim PAH features: 



SFR [M /yr] = 1.18 x 1(T 41 L PAH [erg a" 1 ]. (3) 

A positive correlation between the PAH and neon emis- 
sion is also observed for CSOs (Figure [12]). While a 
slope fitted to the data is consistent with a linear re- 
lation (1.0 ± 0.4), the scatter is nearly an order of mag- 
nitude larger than that measured for the ULIRG sam- 
ple (1.02 ± 0.05). This may be an indication that some 
CSOs have neon emission that does not arise from recent 



star formation. Deviations from the linea r relation also 
are expected for low-met allicity galaxies (jPeeters et al.l 
120041: iCalzetti et ail 120071 ) for which both neon and/or 
PAH features can be significantly depleted with respect 
to the true SFR. 4C +12.50, in particular, deviates sig- 
nificantly both from the ULIRG fit and the remainder of 
the CSO sample (showing almost an order of magnitude 
less PAH emission than would be expected from the neon 
strengths) . 

Star formation rates in CSOs vary by more than two 
orders of magnitude over the sample (Table [7]) , with the 
lowest measured rate in 1146+59 (0.3 — 0.5M©/yr) and 
the highest in OQ 208 and 4C +12.50 (< 45M Q /yr). 
The SFR measured from the PAH strengths is consis- 
tently lower than the SFR from neon by a factor of 4-5; 
however, the conversion between neon and PAH emis- 
sion depends on a number of assumptions (including the 
timescale of star formation and the conversion to Ly- 
man continuum flux). In addition, the data used to cal- 
ibrate Equation [3] is based on data from ULIRGs, with 
unclassified contributions to their active phases; this may 
contribute to the scatter in Figure [12] While emission 
from AGN-heated gas can contribute to neon emission 
(which may explain the large disparity of the SFRs for 
4C +12.50), PAHs are easily dissociated by hard radia- 
tion environments and are almost never detected in AGN. 
For this reason, we regard PAHs as the more robust SFR 
indicator; this also implies that the majority of the star 
formation being detected takes place in the outer regions 
of the CSOs. 

5.5. Black hole masses 

High-ionization transitions in the mid-IR that are as- 
sumed to emit from narrow-line regions gravitationally 
bound to an AGN can b e used to estimate the mass 
of the central black hole (jDasvra et all 120081 ). For the 
six CSOs in which the [O IV] 25.8 /j,m transition is de- 
tected, the line is spectrally resolved for two objects: 
1146+59 and PKS 1718-649. The [O TV]-M BH cali- 
bration yields black hole masses of \og(MsH /Mq) = 9.2 
and 8.6. 4C +31.04 has a measured [O IV] dispersion 
above the instrument al resolution but no t satisfying the 
resolution criterion of lDasvra et al.l ([2008); if valid, it has 
a SMBH mass of log(M Bff /M Q ) = 8.2. The velocity dis- 
persions of two other CSOs in which [O IV] was detected, 
OQ 208 and 4C +37.11, are below the resolution of the 
HR module and give only upper limits on a central black 
hole mass (Table [7]). 

The [O IV] line is also detected in 4C +12.50 with a 
FWHM of 1026 ± 97 km/s, which would correspond to 
a central BH mass of \og(MsH /M Q ) = 11.5 if caused 
entirely by the black hole. This is more than two orders 
of magnitude larger than t he maximum BH ma ss used 
to calibrate the relation of Das vra et al.l ()2008| ) and is 
significantly higher than predicted upper limits on ultra- 
massive black holes. 4C +12.50 has two bright optical 
nuclei within a common envelope separate d by only 1.8" 
(jGilmore fc Shawl llQSfit lAxon et al.N2000T ). which is sev- 
eral times smaller than the width of the LH slit. The 
broad linewidth is likely due to emission from both nu- 
clei moving at relative velocities of a few hundred km/s. 
Several other lines in this object, including [S IV] and 
[Ne III], also have FWHM in excess of 1000 km/s. The 
[O IV] line in this CSO thus cannot be used to estimate 
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Fig. 13.— Distribution of the [Ne IIl]/[Ne II] line ratio; higher 
ratios trace a harder ionizing environment. For comparison, we also 
show neon ratios for samples o f known ULIRGs (Farrah et al. 2007; 
Wille tt et al. 2010), AGN USturm et al.l [2002: JTommasin et al.l 
1 200SI1 . and starburst galaxies llVerma et al.M2003lT ! Dashed lines 
show the mean value for each class of objects. 

an accurate central BH mass. 

Since [O IV] can also be produced in extreme star- 
bursts, the measured line dispersions could represent mo- 
tions of the gas in the ionized environments surrounding 
young, hot stars or the relative motions between star 
forming regions. The star formation rates from 1146+59 
and PKS 1718—649, however, are the lowest measured 
in the sample (< 2 Mq/jt, from both neon and PAH 
data). If star formation is not the main contributor to 
[O IV] emission, then the lines are more likely caused by 
emission from a narrow-line region. Uncertainties in the 
physical location of the ionized gas with respect to the 
black hole, however, contributes to the uncertainty. 

We use a second method to estimate the black hole 
masses from the well-known bulge-luminosity relation. 
We assume that the total luminosity of the galaxy comes 
from the bulge (which could be an overestimate for spi- 
ral hosts or mergers with extended features). We use g 
and r photometry from the SDSS for five CSOs to cal- 
culate Ly', for the remaining galaxies not in the SDSS, 
we use H-band NICMOS and V-band POSS measure- 
ments and apply typical colors for an elliptical galaxy. 
We use the AGN - calibr ated Mbh — Lb u i ge relationship 
from lBentz et al.l ([2009) to compute central BH masses; 
masses for all eight CSOs lie between log Mbh = 8.0 — 
9.0. The measurement for PKS 1718—649 agrees well 
with the mass from the [O IV] dispersion, but the bulge 
Mbh for 1146+59 is lower than the [O IV] estimate by 
~ 1 dex. 

5.6. Ionization and excitation states 



The hardness of the ionizing radiation field in the 300- 
600 A range can also be traced using the excitation states 
of fine-structure lines from the same species; a harder 
radiation field increases the line ratio of higher/lower 
states. In the mid-IR, the [Ne III]/[Ne II] ratio is the 
only pair of lines visible in a majority of the CSOs. We 
plot the line ratios for the eight objects in which neon is 
detected (Figure [T3|) . along with the distribution of other 
IR-bright galaxies measured with the ISO and Spitzer. 

The neon ratios in all CSOs are between 0.4 and 
1.0 . The aver a ge neo n rat io from the AGN samp les 
of iSturm et all (|2002Tl and iTommasin et all (120081) is 
1.3 + _ 1.0. Starbursts with mid-IR data (jVerma et all 
2003) are much fewer in number and show a larger spread 
in neon ratios, with a mean value of 1.8 + 3.3. 

The distribution of CSO neon ratios is not significantly 
different from either starbursts or AGN, although com- 
parisons are limited due to the small number of CSOs. 
The mean hardness ratio in AGN is higher than that 
in the wider ULIRG population, likely due to the harder 
X-ray/UV environment that can arise from accreting gas 
around a central black hole. Given the large overlap in 
hardness ratios for both AGN and starburst galaxies, this 
parameter alone cannot classify CSOs in either category. 

5.7. Mid-IR classification of the CSO hosts 

The mid-IR spectra show a heterogeneous set of fea- 
tures for the CSOs, as discussed above. While the objects 
were selected based on similar characteristics in the ra- 
dio, analysis of the mid-IR features shows significant de- 
viations from a single CSO "family" . We classify the ob- 
jects into four categories, based on features visible in the 
LR spectra (PAH strength, silicate depths, fine-structure 
and H2 lines): 

• Category 1 shows several weak-to-moderate PAH 
emission features, often with only an upper limit 
on the 6.2 fim. feature, with moderate 6*9.7 absorp- 
tion. This is consistent with low levels of star for- 
mation and a moderately dusty environment. This 
is the most common category for the CSOs ob- 
served, including 4C +31.04, 1146+59, 4C +12.50, 
PKS 1718-649, and 1946+70. 

• Category 2 is similar to objects in Category 1, with 
weak PAH emission and moderate S9.7 absorp- 
tion; however, these objects display strong H 2 fea- 
tures that may arise from shocked molecular gas. 
4C +37.11 is the lone example of this in the sample. 

• Category 3 also has weak PAH emission on top 
of a strong continuum, but shows the silicate dust 
features in emission rather than absorption. This is 
a strong signatur e of a quasar or AGN in the galaxy 
(jHao et al.l 120071 ) . although active star formation 
may still be taking place. OQ 208 is the only object 
in this category. 

• Category 4 objects have no discernible emission 
features, with the only break in a continuum that 
can be fit with a single power law the silicate ab- 
sorption at 9.7 and 18 fxm. PKS 1413+135 is the 
only CSO in the sample with these basic blazar 
characteristics. 
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The IRS non-detection of the CSO 1245+676 (§ ED} 
is also interesting. The object is detected with 2MASS 
with rriK, = 13.2, suggesting the existence of a large 
population of old, red stars. The low levels of mid-IR 
emission, however, imply that the object is dust-poor 
and has no strong PAH features which would signal high 
levels of star formation. These observations are consis- 
tent with the standard picture of an AGN host as an 
elliptical galaxy with little gas or dust. The host galaxy 
of 1245+676 thus might be an example of a dust-poor 
galaxy similar to Categories 1 or 2; deeper observations 
that can place strong limits on the PAH emission and sil- 
icate absorption are critical to answering this question. 

In summary, 5/8 CSOs have moderately dusty envi- 
ronments and moderate levels of star formation. The 
exceptions from this typically show strong evidence for 
being AGN-dominated, including the broad-line galaxy 
OQ 208 and the BL Lac PKS 1413+135. 

6. DISCUSSION 
6.1. Photoionization vs. star formation 

One of the key issues in interpreting the mid-IR data 
is the ionization source of the emission features. For IR- 
bright galaxies, this typically has two main possibilites: 
photoionization of the gas by an AGN, or from star for- 
mation possibly spread throughout the galaxy. If the 
former is the case, then it is likely that some of the gas 
is close enough to accrete onto the AGN and provide a 
power source for the radio jets in the inner kpc. 

For several CSOs, hard X-ray emission provides direct 
evidence for supermassive BHs; in the case of 4C +37.11, 
a pair of BHs has been directly detected. In the op- 
tical, complex Balmer profiles in OQ 208 also identify 
a broad-line region surrounding an AGN. The fork dia- 
gram places all eight CSOs on the AGN branch, with the 
silicate absorption ratios favoring a clumpy dust torus for 
most galaxies in the sample. Finally, multiple CSOs show 
very low levels of star formation (from multiple tracers) . 
PKS 1413+135 has only upper limits on star formation 
from mid-IR tracers, although the limits are higher than 
most SFR measurements in the sample due to its much 
greater distance. 

There is also evidence that star formation may be re- 
sponsible for a significant fraction of the mid-IR emis- 
sion, however. Two CSOs (4C +12.50 and OQ 208) 
show strong PAH features that would normally be dis- 
sociated if directly exposed to an AGN, with Lju and 
Ljve also supporting SFR > 50M Q /yr in both galax- 
ies. The detection of H2 in most galaxies indicates the 
presence of molecular gas - this has also been detected at 
cooler temperatures at mm- wavelengths in several CSOs, 
meaning that star-forming gas is still present in most of 
the sample. [Ne V], a "smoking gun" of AGN activity 
in the mid-IR, is not detected in any CSO, while the 
[O IV] emission could arise from photoionization from 
either O and B stars or an AGN. Finally, the poor fits 
of many CSOs to the classic AGN torus of clumpy dust 
could indicate a different geometry, with distributed dust 
and gas undergoing star formation in extended features. 

The results of the mid-IR data also have intriguing 
implications for the evolutionary history of radio galax- 
ies. Although the sample size is small, nearly all CSOs 
show contributions from both AGN and star formation 



(Figs. I5HTU]). This is strong support for the scenario in 
which young radio galaxies are formed by the merger of 
two disks; remnant star formation from the initial merger 
is still detected, but clear signatures of AGN have also 
begun to appear in other galaxies ([O IV], hot dust, and 
a clumpy dust geometry in the mid-IR; X-ray emission 
and broad- line regions also support this). Since we are 
deliberately selecting only the youngest radio galaxies, 
it is not surprising that the AGN do not yet dominate 
the IR emission if they are also in their earliest stages. 
As the galaxies continue to virialize (and possibly feed 
the central BH), we predict that star formation tracers 
such as PAH emission should become weaker, the AGN 
and hard radiation environment will increase, and the 
jets continue t o advance, event ually forming FR II-type 
radio galaxies (Bcgclmanl ll996| ). Mid-IR observations of 
later-stage radio galaxies ("medium symmetric objects", 
or MSOs) could also support this hypothesis if they show 
a larger AGN fraction than their younger counterparts. 

6.2. What is the power source for the radio jets? 

Mid-IR data can discriminate between two quite dif- 
ferent models for powering radio-loud AGN: accre- 
tion from a disk and/or torus as opposed to the tap- 
ping of spin energy from a supermassive black hole 
(|Blandford fc Znaieldll977t IWilson fc ColT^ll995l) . If 
accretion is the energy source, then there should be ev- 
idence for enough dust and gas in the nuclear region to 
obscure the AGN in near-IR NICMOS images without 
obvious broad lines. Furthermore, we would expect to 
detect highly-ionized fine -structure transitions such as 
[Ne V] (|Armus et al.ll2004T) . The slow Bondi accretion vs. 
ADAF models, however, are not so easily distinguished 
from mid-IR data. 

Nearly 75% of all QSOs show [Ne V] emission and 
> 90% show [O IV] 25.9 /im emis sion indicative of accret - 
ing gas close to the black hole dSchweitzer et al.l [20061. 
The sample of lArmus et al.l (|2004[ ) also detected [Ne V] in 
3/10 ULIRGs. If the absence of [Ne V] is a physical re- 
sult (and not a question of sensitivity), then this requires 
a substantial energy source not powered by accretion. 
The non-detection of [Ne V] in all eight CSOs, based on 
simple Gaussian statistics for QSOs, is significant at the 
~ 4<7 level. 

The lower ionization potential transition [O IV] was 
detected in six CSOs; this line can be associated both 
with AGN and photoionization from extreme starbursts. 
Possible scenarios in which an AGN could produce 
[O IV] without [Ne V] are: significant depletion of neon 
in the accreting g differential extinction screen that 
obscures the neon emission; or, a radiation environment 
in which the majority of the ionizing flux has energies be- 
tween 55 < Ei on < 97 eV. This may be the case in lower- 
luminosity AGN, compared to luminous QSOs. The lim- 
its on [O IV] and [Ne V] 24 /*m for the most continuum- 
dominated CSOs (PKS 1413+135 and 4C +12.50) are 
nearly an order of magnitude higher than the bright- 
est detections. We also test the [Ne V] upper limits by 
comparing it to the average [Ne V]/[Ne II] ratios in PG 
quasars, which are robust for 6/7 CSOs with [Ne II] de- 
tections and indicating that signal-to-noise ratio is not 
the primary factor in the absence of [Ne V] . 

We can also address the possibility of accretion fuel- 
ing the radio jets from an energy balance argument. The 
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power in the radio jets has two components: the mechan- 
ical power necessary to excavate the radio cavity in the 
ISM and the radio luminosi ty from the jets themselves. 
Following I Allen et al.l (|2006| ). the work needed to expand 
the cavity is: 

Pjet = -^PV/tage, (4) 

where p = kTp/ p,m p is the pressure of the surround- 
ing ISM and 7 is the adiabatic index of the gas (as- 
sumed to be 5/3). We assume a temperature of 1 keV 
and n ~ 1 cm" 3 , based on the largely homo geneous re- 
sults f or elliptical galaxies in the sample of lAllen et al.l 
(2006). While the cavities cannot be directly observed 
without high spatial resolution X-ray images, we use the 
projected dimensions of the radio jets as an estimate. We 
assume V — Ttrir^ is the volume of a cylindrical cavity 
with projected radii of r w and n along the jet minor and 
major axes. We combine r; with the kinematic ages (t age ) 
of the CSOs from VLBI proper motions to determine the 
timescale necessary to excavate the cavity. 

We add Pjet to the total radio luminosity assum- 
ing a self-absorbed synchrotron spectrum with the low- 
frequency emission S oc v^l 2 . The high-frequency spec- 
tral index is measured (where possible) using the 4- 
8 GHz fluxes. In almost all CSOs, we find that P ra dio is 
2-3 orders of magnitude larger than Pjet , the opposite of 
jet-dominated FR I galaxies. 

To test the accretion model, we calculate the Bondi 
power for spherically symmetric absorption: 

PBondi = 4TT\G 2 Ml H pr 1C 2 c- 3 , (5) 

where A = 0.25 is a numerical coefficient depending on 
the adiabatic index of the accreting gas, 77 is the effi- 
ciency of accretion, c s and p are the sound speed and 
density of the gas, and Mbh the mass of the centra l 
black hole onto which the gas is accreting (Bondi 1952). 
The number density i s consistent with the findings of 
iPihlstrom et al.l (|2003f ) requiring Uism ~ 1 — 10 cm -3 
to confine GPS/CSS sources. The black hole mass is 
estimated using the bulge luminosity- Mbh relation. 

Models o f jet power both in ellipt ical galaxies 
(I Allen et al.l [2006) and galaxy clusters (|Birzan et al.l 
[200¥ demonstrate a scaling relation between Bondi ac- 
cretion and the jet power. Results for CSOs are shown 
in Figure [HI For an assumed accretion efficiency of 
r\ = 0.1, Bondi accretion is capable of powering the radio 
jets in most objects, with the prominent exceptions of 
4C +12.50 and PKS 1413+135. If advection-dominated 
accretion flows can be shown to exist without produc- 
tion of the high-ionization [Ne V] lines, then we must 
still consider gas accretion as a viable power source. 

As an alternative, we examine the evidence for the 
iWilson fc Colbertl ([1995ft model in which the spin energy 
of two SMBHs is tapped to power the radio jets in the 
absence of accreting gas. Since the delay between the 
merger and the onset of radio activity and the timescale 
of merging black holes is both ~ 10 8 yr. However, the 
general model offers relatively few testable predictions 
that can be made in the absence of accreting gas. One 
constraint, however, is that the BHs must be of compa- 
rable size and have masses greater than 10 8 Af Q . In the 
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Fig. 14. — Bondi accretion power as a function of the combined 
jet mechanical power and radio luminosity. The red dots are CSOs, 
while the b l ack d ots are X-ray luminous elliptical galaxies from 
lAllen et al.l Q2006f ). Arrows show the limits on CSOs without con- 
strained tage- The dashed line is a linear fit to the X-ray ellipticals; 
we label CSOs that deviate significantly from this fit. 

case of 4C +37.11, these conditions have already been 
met; the BH masses measured both with Lb u i ge and the 
[O IV] dispersion also lie above this line. The separation 
of the two BHs (7.3 pc), however, mean that the cores 
are still too distant to tap the spin energy for powering 
the radio luminosity. No radio-loud galaxy has yet been 
discovered in which merging BHs are close enough to tap 
their combined spin energy. 

7. CONCLUSIONS 

We observed a sample of eight \aw-z CSOs with the 
IRS, representing very young radio galaxies and possi- 
ble progenitors to FR lis. Spitzer investigations have 
revealed: 

1. The mid-IR SEDs for CSOs show significant va- 
riety. This includes silicate features seen in both 
emission and absorption, spectral indices indicat- 
ing dust at a wide range of temperatures, and fine- 
structure and PAH emission. We sort the galax- 
ies into four categories based on their IRS spec- 
tra, with the most common features including weak 
9.7 /xm dust absorption, moderate PAH and atomic 
emission, and H2 temperatures of 200-400 K. There 
is no evidence in the mid-IR, however, of a single 
"family" for CSOs. 

2. PKS 1413+135 is a unique object among CSOs, 
showing no emission or absorption features associ- 
ated with gas in the ISM. The only identifiable fea- 
tures in the mid-IR are the silicate dust troughs at 
9.7 and 18 /xm. While the dust redshifts are roughly 
consistent with the near-IR obscuration, it does not 
address the possibility that the AGN lies along the 
line of sight and may not be physically associate d 
with the optical counterpart ([Perlman et alj |2002) . 

3. Numerical models indicate that a clumpy dust 
torus can fit the silicate features in most CSO spec- 
tra. Features common to most models have 1-10 
dust clouds along the line of sight, with optical 
depths reaching as high as a few hundred in ry. 
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Multi- wavelength evidence (including X-ray fluxes, 
silicate/PAH ratios, and radio- imaging of SMBHs) 
indicate that significant fractions of the sample 
have AGN at their centers. Evidence for contin- 
uing star formation, however, suggests that CSOs 
occupy a continuum between starburst galaxies and 
AGN. This supports the scenario that the elliptical 
hosts were formed by a merger of disk galaxies with 
a delay before the activation of an AGN. 

Based on the work necessary to expand the radio 
jet cavities, accretion is a viable power source for 
6/8 CSOs in the sample. The non-detection of 
high-ionization lines and absence of large extinc- 
tion support an ADAF model, if this is the case. 
Although BHs with M > 10 8 A/ Q likely exist in all 
CSOs, there is no direct evidence that BH spin en- 



ergy provides the power source for the radio jets. 



This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory (JPL) and the California In- 
stitute of Technology. We also made use of the 
NASA/IPAC Extragalactic Database (NED) which is op- 
erated by JPL and Caltech under contract with NASA. 
Many thanks are due to L. Armus, V. Charmandaris, and 
H. Spoon for their help with data reduction, K. Dasyra 
and B. Keeney for useful discussions, and to the Spitzer 
Science Center for hosting KW and JD while collaborat- 
ing on data analysis. KW and JS gratefully acknowledge 
support from Spitzer programs 30515 and 50591. 



APPENDIX 

In this section, we summarize the multi-wavelength properties of the CSOs from the literature. Objects are referred 
to primarily by their radio nomenclature, rather than the host galaxy. 

4C +31.04: The radio emission has a compact core with symmetric hotspots and an unusual "hole" in its eastern lobe. 
iGiroletti et all (|2003h measure lobe advance speeds of (0.33 ± 0.06)/i 65 c and a kinematic age of 550 ± 100 years. The 
distance from the core to the hotspot edges is 50-60 pc for both lobes. The CSO's host galaxy is MCG 5-4-18, which 
has a nearly edge-on circumnuclear disk and strong ob scuration in both the nucleus and outer regions. "Cone-like" 
features in the nucleus, seen bot h in H i (jConwavHiggl ) and HST images (jPerlman et al.ll2001bl ). are aligned with the 
radio jets. iGarcia-Burillo et al.1 (j2007j ) detect a molecular /dusty disk around the AGN with M gas > 5 x 10 9 M© and 
a size of 1.4 kpc. Distortions and non-circular motions in the molecular disk suggest it is not dynamically relaxed, 
possibly resulting from either a recent merger, gas accretion events, or interactions between the radio jet and the disk. 
4C +37.11: This galaxy possesses several features highly unusual among CSOs. It has extended kiloparsec radio 
emission that may be either due to a merger or previous activity in one of the cores. The parsec-scale radio structure 
shows two compact, flat-spectrum sources identified as supermassive black holes (SMBHs) that have not yet coalesced 
(M aness et "aLll2004f ). with a combined mass of 1.5 x 10 8 M© and a separation of only 7.3 pc dRodriguez et al.l | 2006[) . 
The kinematic signature of the BH also appears in a pair of broad Hi absorption features ()Morganti et al.l 2009) . 



The radio jets are 15-25 pc in length, with the north hotspot moving outwards at ( 0.137 ± 0.034) / i7 1 1 c; this yields a 
kinematic age of 502 ± 129 years. The host galaxy is the elliptical 0402+379, which IStickel et all (|1993f) describe as 
having a "flat brightness distribution" with no central point source in the optical. The X-ray brightness of this galaxy 
is the highest in the sample, and it is the only known CSO detected in the ROSAT All-Sky Survey. 
1146+59: This is the nearest and least radio-luminous galaxy in the sample. It has a double-lobed radio continuum 
stru cture with advance sp eeds of (0.23 ± 0.05)/ij" O oC (|Tavlor et al.lll998l ). with Hi detected in absorption against both 
jets (|Peck fc T aylor 1998]). The separation speed between the jets is (0.23 ± 0.05)c with a separation of only ~ 2 pc; 
if constant, this means the jets emerged only 34 years ago, making it the youngest known CSO. The host galaxy, 
NGC 3894, is classified as an E/S0 with an optical absolute magnitude of approximately L*. HST observations with 
WFPC2 reveal two du st lanes nearly perpen dicular to the radio axes, interpreted as the halves of a dusty torus at a 
high inclination angle (jPerlman et al.ll2001bO . The galaxy was observed with Chandra in early 2009 (Perlman et al., 
in prep). 

4C +12.50: The northern radio jet in this galaxy is much fainter and shorter than the southern counterpart, which 
extends out to ~ 150 pc. There is a s mall amount of kpc-scale ra dio emission suggesting that the galaxy was active 
before the current burst of jet activity ([Condon fc Brodericklll99"fl) . The velocities of the northern jet components are 
(1.0 ± 0.3)/i 7 Tq 1 c and (1.2 ± Q.2)h^c, representing one of the only known examples of superluminal motion in CSOs; 
the velocity of the sout hern hotspot, however, is poorly constrained. Assuming a "typical" CSO advance speed of 
0.3c. iLister et al.l ()2003[ ) calculate a kinematic age of 1700 yr; th is is consistent with the dynamics of a precessing jet, 
requiring an age less than 10 5 yr. iDarling fc Giovanellil <|2002f ) also report a tentative detection of OH megamaser 
emission. Optical images of the host galaxy, IRAS 13451+1232, show irregular isophotes in a double nucleus, a 
tidal tail, and several companion galaxies indica t ing a n advanced merger with little star formation (|Gilmore fc Shawl 
119861: ISurace et aLlll998l:ISurace fc Sandersl[l999l l200"ol ). The low levels of star formation are puzz ling, given the large 
amounts of molec ular gas (6.5 x 10 10 M©) detected near the nucleus and its high IR luminosity (iMirabel et al.l[l989t 
lEvans et al.lll999h . High amounts o f reddening and the presence of hard X-rays are interpreted as a quasar obscured 
behind a foreground dust screen bv lO'Dea et all (l2000Tl . 

OQ 208: VLBI observations show a slightly asymmetric distribut ion between two radio l obes at a maximum sep- 
aration of ~ 8 pc, with no strong emission visible from the core (Stangh ellini et aTlll997| ). The spectral turnover 
can be explained both by free-free and synchrotron self- absorption models, with lobe advance speeds of 0.33/iyo e 
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(|Kameno et al.l [20001 : Kie et al.l l2~005) and a kinem atic age of 92 yr. The h ost galaxy, Mrk 668, shows the presence of 
both a tail and galactic companions in the optical (Stangh ellini et aLlll993l ). suggesting that the galaxy is dynamically 
distu r bed. Complex Ba l mer li ne profiles classify the host galaxy as a broad line radio galaxy, or BLRG ()Marziani et al.l 
f993). iGuainazzi et al.l (|2004[ ) detect flat X-ray emission from 2-10 keV with powerful K a iron line emission, the first 
known example of a Compton-thick AGN in a BLRG. This is the only CSO in the sample with clear evidence of a 
broad-line region. 

PKS 1413+135: This is the most radio-luminous and by far the most distant galaxy in the sample. At milliarcsecond 
scales, the galaxy has both a jet and counterjet with a strong core component, with hotspots extending out to ~ 139 pc 
for the counterjet. The jets are likely oriente d close to the line of sight ; pressure balance arg uments suggest an age 
of 10 3 - 10 4 yr and velocities of ~ 0.05/i^c (jPerlman et al.lll996l 120021 ) . iStocke et all (|1992| ) suggest that the radio 
component may be background to the optical counterpart, based on the absence of reprocessed radiation from near-IR 
dust continuum and emission lines. The host is a spiral galaxy wit h prominent dust lan es , a nuclear point source in 
the NIR, and a blue luminosity consistent with an Sbc classification (|Perlman et al. 2002). iWiklind fc Combei ()1997f ) 
detect several absorption lines in the spiral galaxy at z = 0.247, consistent wit h a two- phase medium with dense 
molecular clumps embedded in a more diffuse component. iPerlman et al.l (|2002[ ) suggest that the absorbing screen 
consists of a giant mole cular cloud (G MC) within the outer disk of the optical galaxy; this is consistent with OH 
absorption measured by I Darling ()2004l ). 

PKS 1718—649: The radio emission has two lobes separated by ~ 2 pc and no apparent emission from a core 
(jTingav et alj [2002). The upper limit on the lobe separation speed is 0. 08/i 100 c; iTingav et al.l (|1997f ) adopt a kinematic 
age limit of < 10 5 yr based on the jet energy and density advancing into an ISM with n < 3 x 10 6 cm -3 . The majority 
of the Hi is distributed in an incomplete ring around the nucleus with a diameter of 37 kpc; ki nematics of a weak 
Hi envelope extending out to 180 kpc suggest a past merger with at least one gas-rich spiral galaxy ()Veron-Cettv et ahl 
119951) . The total Hi mass from VLBI is 3.1 x 10 10 M W . T he host galaxy is NGC 6328, a LINER with photoionization 
as the dominant excitation mec hanism (lFilippenko]J1985l) . Optical images show a bright elliptical nuclear region and 
faint, extended spiral-like arms (Vcron- Cettv et al.lll995D . A bar-like structure in the nuclear emission region imaged 
in Ha and [N II] is coincident with the compact radio source ([Keel fc Windhorstl[l99ll) . 

1946+70: The milliarcsecond radio st ructure shows an inverted core and two jets with multiple hotspots and "S"- 
symmetry extending out to 25-35 pc (|Tavlor fc Vermeulenlll997[). The expans ion velocity between the hotspots is 
0.024c, corresponding to a kinematic age of 4000 ± 1000 years ( Taylor et al.ll200"9l ). Ejection of pairs of jet components 
appear to take place on timescales of ~ 10 yr, with the most recent outburst observed in 1997. The peak Hi absorption 
occurs nea r the c ore and extends outward somewhat into the jets, suggesting the presence of a foreground screen. 
IPeck et al~l (|1999D propose a model in which the central jets are obscured by a torus of clumpy gas (both neutral and 
ionized) at 50-100 pc, with a warped molecular disk surrounding the AGN at ~ 1 pc. HST observations of the ho st 
galaxy show mild isophotal twists and ellipticities, which may be indicative of a past merger (Pcrlma n et al.ll2001bD . 
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